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We  use  a  two-dimensional  nonlinear  compressible  mesoscale  model  to  study 
mountain  waves  in  the  middle  atmosphere  over  the  southern  Andes.  Using  re¬ 
alistic  temperatures,  winds  and  topography,  the  model  generates  large-amplitude 
long-wavelength  breaking  mountain  waves  in  the  middle  atmosphere  that  compare 
favorably  with  satellite  measurements.  Modeled  waves  have  preferred  horizontal 
wavelengths.  Spectral  analysis  reveals  correspondences  between  these  wavelengths 
and  peaks  in  the  spectrum  of  Andean  topographic  elevations.  The  shorter  waves 
reach  the  stratosphere  well  before  the  longer  ones,  consistent  with  group  veloc¬ 
ity  arguments,  with  longer  wavelengths  ultimately  dominating.  At  later  times  we 
find  evidence  of  downward  propagating  secondary  waves  produced  by  upper-level 


breaking  of  the  primary  waves. 

1.  INTRODUCTION 

Theory  and  parameterization  have  indicated  that  gravity 
waves  forced  by  flow  over  mountains  are  a  major  source 
of  dynamical  driving  for  the  stratosphere  and  mesosphere 
[e.g.,  Bacmeister,  1993;  Boville,  1995).  Detailed  mesoscale 
model  simulations  have  proved  useful  in  modeling  mountain 
wave  dynamics  and  refining  simpler  parameterizations  of 
these  processes  in  the  troposphere  and  lower  stratosphere 
[e.g.,  Lap  rise,  1993;  Kim  and  Arakawa,  1995;  Durran,  1995; 
Broad,  1996;  Do  mb  rack  et  al„  1998]. 

However,  only  a  few  such  models  have  incorporated  a 
middle  atmosphere  and  investigated  mountain  wave  dynam¬ 
ics  up  to  and  beyond  the  stratopause.  One  of  the  first  rela¬ 
tively  detailed  studies  was  by  Schoeberl  [1985],  who  derived 
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linear  nonhydrostatic  mountain  wave  solutions  within  real¬ 
istic  wind  profiles  from  the  ground  up  to  ~70-80  km.  More 
detailed  numerical  simulations  were  performed  by  Bacmeis¬ 
ter  and  Schoeberl  [1989]  using  a  two-dimensional  nonlinear 
time-dependent  anelastic  model  of  flow  over  an  obstacle. 
They  found  zones  of  strong  wave  overturning  in  the  middle 
atmosphere  that  generated  secondary  downward-propagating 
waves.  Their  model  used  idealized  representations  of  the 
background  atmosphere.  Satomura  and  Sato  [1999]  used  a 
fully  nonlinear  compressible  two-dimensional  model  of  flow 
over  orography  with  realistic  representations  of  the  lower 
and  middle  atmospheres,  based  on  climatology.  They  too 
simulated  strong  wave  breaking  in  the  middle  stratosphere 
and  generation  of  secondary  waves  from  these  unstable  re¬ 
gions,  although  the  wavelengths  differed  from  those  found  by 
Bacmeister  and  Schoeberl  [1989].  All  of  these  studies  con¬ 
sidered  flow  over  a  bell-shaped  obstacle  and  omitted  Coriolis 
accelerations. 

Here,  we  take  further  steps  towards  “real  world"  mesoscale 
model  simulations  of  mountain  waves  in  the  middle  atmo¬ 
sphere.  Our  simulations  focus  on  a  specific  set  of  satellite 
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measurements  of  mountain  waves  over  the  region  shown  in 
Plate  la.  Colored  squares,  labeled  1,2  and  3,  show  locations 
of  three  successive  vertical  temperature  profiles  acquired  by 
the  Cryogenic  Infrared  Spectrometers  and  Telescopes  for  the 
Atmosphere  (CRISTA)  on  6  November,  1994  at  6:24  UT. 
These  profiles,  reproduced  in  Plate  4  of  Grossman  [2000], 
reveal  coherent  stratospheric  temperature  oscillations.  Eck- 
ermann  and  Preusse  [1999]  used  theory  and  parameteriza¬ 
tion  to  argue  that  these  oscillations  were  large-scale  mountain 
waves  generated  by  flow  over  the  Andes. 

In  section  2,  we  show  that  two-dimensional  mesoscale 
model  simulations  offer  a  potentially  useful  means  of  study¬ 
ing  the  generation  and  evolution  of  the  mountain  waves  ob¬ 
served  on  this  day.  We  describe  our  mesoscale  model  in 
section  3,  and  apply  it  in  section  4  using  “best"  estimates 
of  the  large-scale  winds,  temperatures  and  orography.  The 
results  are  analyzed  in  section  5  and  the  findings  are  summa¬ 
rized  in  section  6. 

2.  CRISTA  OBSERVATIONS  AND  REGIONAL 
DOMAIN 

Plate  1  shows  that  the  large-scale  topography  here  is  quasi- 
two  dimensional.  The  “ridge  axis”  in  Plate  la  is  a  subjective 
fit  to  the  orientation  of  the  long  axis  of  the  Andes.  The  “per¬ 
pendicular  axis,”  aligned  normal  to  this  ridge  axis,  passes  as 
close  as  possible  to  the  three  CRISTA  profiles,  and  illustrates 
that  the  CRISTA  data  were  acquired  along  an  orbital  segment 
that  passed  nearly  orthogonally  across  the  ridge  axis  of  the 
Andes. 

We  take  the  blue  dot  in  Plate  1  at  the  western  end  of  the 
perpendicular  axis  as  our  “upstream  point,"  since  surface 
winds  were  eastward  and  the  atmosphere  here  is  far  enough 
upstream  to  remain  unaffected  by  the  Andes  [Seluchi  et  al.. 
1998].  Figure  la  plots  a  hodograph  of  the  horizontal  wind 
profile  at  the  upstream  point  on  6  November,  1994  at  6:00  UT, 
using  data  from  NASA's  Data  Assimilation  Office  [ Coy  and 
Swinbank,  1997].  Upstream  troposphere-stratosphere  winds 
were  roughly  parallel  to  the  perpendicular  axis  at  this  time, 
particularly  near  the  surface  and  in  the  stratosphere,  our  pri¬ 
mary  regions  of  interest.  Maximum  deviations  of  ~45°  oc¬ 
curred  near  the  tropopause  in  Figure  la. 

Similar  conditions  occur  when  eastward  surface  winds 
flow  in  from  the  ocean  across  the  Norwegian  Mountains. 
In  these  situations,  three-dimensional  mesoscale  models  of¬ 
ten  predict  long  two-dimensional  (plane)  mountain  waves 
in  the  lower  stratosphere  over  Scandinavia,  which  compare 
well  with  aircraft  and  radiosonde  data  [Leutbecher  and  Volk- 
ert ,  1996;  Dombrack  et  al.,  1998,  1999].  Consequently, 
these  waves  can  be  modeled  fairly  accurately  using  a  two- 
dimensional  mesoscale  model  [Volkert and Intes,  1992].  We 
too  use  a  two-dimensional  mesoscale  model  as  a  first  step  in 
this  study. 


3.  NUMERICAL  MODEL  AND  DOMAIN 
3.1.  Model  Description 

The  model  used  here  was  developed  from  an  earl  ier  version 
first  described  by  Tan  and  Leslie  [1998].  It  solves  the  Eu¬ 
ler  equations  for  fully  compressible  nonlinear  nonhydrostatic 
atmospheric  flow  over  arbitrary  terrain  in  a  two-dimensional 
spatial  domain.  The  velocity  u  =  (u,t;,  w),  potential  tem¬ 
perature  0  and  Exner  pressure  ir  are  split  into  background 
and  perturbation  terms,  as  follows: 

«(*.*.<)  =  u(z)  +  u'(ar,  z,  t), 

9{x,z,t)  =  0(z)  +  O'(x,z,t), 

tr(z,z,<)  =  ir(z)  +  x'(x,z,t). 

The  velocity  components  v  and  w  contain  perturbation  terms 
only.  A  background  temperature  profile  T(z)  is  used  to 
initialize  the  background  density  p(z),  potential  temper- 
aure  0(z)  and  Exner  pressure  # (z)  =  T(z)/0(z),  assuming 
hydrostatic  equilibrium.  The  cross-ridge  background  flow, 
ti(z),  is  assumed  to  be  in  thermal  wind  balance.  The  mo¬ 
mentum,  thermodynamic  and  continuity  equations  become: 


du'  ,du  dir'  , 

=  D(u'), 
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cp  and  c„  are  the  mass  specific  heats  at  constant  pressure  and 
volume,  respectively,  R  =  cp  —  cv  is  the  gas  constant  for 
dry  air,  g  is  the  gravitational  acceleration,  c,  is  the  speed  of 
sound  and  /  is  the  Coriolis  parameter.  The  D(X')  terms 
in  (laMld)  are  damping  terms  for  the  perturbation  quanti¬ 
ties  X',  and  represent  the  combined  effects  of  a  first-order 
turbulence  closure  scheme  [Lilly,  1962],  fourth-order  compu¬ 
tational  damping  (to  suppress  small-scale  numerical  modes). 
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cause  numerical  aliasing.  In  all  experiments,  wind  profiles 
were  ramped  up  over  2  hours  from  an  initial  state  of  rest,  to 
minimize  the  generation  of  transient  oscillations  [e.g.,  Ikawa, 
1990].  The  upper-level  sponge  was  effective  above  ~55  km. 
Equations  were  forwarded  using  Ati  =  5  s  and  At,  =  1  s. 

4.  MODEL  SIMULATIONS 

4.1.  Experiments  A,  B  and  C 


Figure  1.  (a)  Hodograph  of  horizontal  winds  over  the  upstream 

point  in  Figure  1 ,  taken  from  DAO  assimilated  data  on  6  November. 
1994,  6:00  UT.  The  altitudes  (in  kilometers)  of  each  wind  value 
arc  labeled  along  the  hodograph  trace.  Orientation  of  the  perpen¬ 
dicular  axis  in  Plate  la  is  shown  with  a  broken  line;  (b)  black 
solid  curve  shows  the  upstream  DAO  winds  along  the  perpendic¬ 
ular  axis,  i(z).  Other  labeled  curves  show  model  u(«)  profiles 
used  in  Experiments  A.  B  and  C;  (c)  vertical  profile  of  upstream 
DAO  temperatures  (solid  black  curve),  and  the  model  curves  used 
in  Experiments  A,  B  and  C. 


and  a  sponge  layer  (Rayleigh  damping)  at  upper  levels  [Dur- 
ran  and  Klemp.  1983]. 

To  incorporate  topography  at  the  bottom  boundary,  the 
model  uses  the  terrain-following  vertical  coordinate 


L[z-h(x)\ 

L-h(x) 


(2) 


where  L  is  the  top  model  altitude  and  h(x)  is  the  topographic 
height  function  [Gal-Chen  and Sommenille,  1975].  A  free- 
slip  lower  boundary  condition  is  used. 

The  model  equations  (la)-(le)  are  solved  numerically 
using  the  Klemp- Wilhelmson  time-splitting  method  on  an 
Arakawa  C-grid  [ Klemp  and  Wilhelmson,  1978;  Durran  and 
Klemp,  1983;  Wicker  and  Skamanock,  1998].  The  second- 
order  Runge-Kutta  method  is  used  for  the  large  time  steps, 
Ati,  and  the  forward-backward  scheme  is  used  for  the  small 
time  steps,  Al,.  An  implicit  method  is  used  in  the  verti¬ 
cal  for  the  acoustic  terms  in  the  governing  equations.  A 
third-order  upwind  finite  difference  scheme  is  used  for  the 
advection  terms  and  second-order  centered  finite  differenc¬ 
ing  is  used  for  the  pressure  gradient  and  divergence  terms. 
The  Miller-Thorpe  radiation  condition  is  used  at  the  lateral 
boundaries  to  prevent  side  reflections  [ Miller  and  Thorpe, 
1981].  The  current  model  configuration  was  validated  by 
reproducing  mountain  wave  patterns  from  previous  two- 
dimensional  model  experiments  [e.g.,  Doyle  et  al.,  2000]. 

3.2.  Model  Domain  and  Initialization 

The  perpendicular  axis  x  in  Plate  1  extended  to  3000  km 
in  our  model  with  a  grid  resolution  Ax  =  2  km.  We  set  the 
model  lid  L  =  80  km  and  vertical  resolution  Az  =  0.5  km. 
Topographic  elevations  h(x)  in  Plate  lb  were  smoothed  to 
remove  the  two  and  four  grid-point  wavelengths  that  can 


Three  different  model  experiments,  denoted  A,  B  and  C, 
were  conducted.  As  shown  in  Figures  lb  and  lc,  each 
successive  experiment  more  closely  approximated  the  up¬ 
stream  DAO  wind  and  temperature  profiles  at  the  time  of  the 
CRISTA  measurements. 

Experiment  A  used  an  idealized  wind  profile  u(z)  that 
omitted  the  westward  shear  layer  in  the  DAO  profile  above 
35  km  (Figure  lb).  This  control  experiment  allowed  us  to 
study  mountain  wave  generation  and  propagation  without  the 
complication  of  critical-level  effects  (u(x)  =  0). 

In  Experiment  B,  a  linear  westward  shear  layer  was  intro¬ 
duced  above  35  km,  yielding  ti(x)  =  0  at  z  «  43  km.  Experi¬ 
ment  B  allowed  us  to  introduce  the  critical  layer,  while  main¬ 
taining  an  identical  atmospheric  situation  to  Experiment  A 
below  35  km.  To  this  end,  both  experiments  used  background 
temperatures  'T(x)  from  the  1976  U.S.  Standard  Atmosphere 
at  mid-latitudes  [ Minvier ,  1977],  which  fit  the  upstream  DAO 
temperatures  from  ~0-55  km  quite  well  (Figure  lc). 

Experiment  C  used  cubic  spline  fits  to  the  DAO  wind  and 
temperature  profiles  from  0-55  km.  This  yielded  a  zero  wind 
line  at  z  k  46  km,  slightly  higher  than  in  Experiment  B. 
Lacking  DAO  data  above  55  km,  from  55-80  km  we  used 
constant  winds  and  a  linear  temperature  gradient  equal  to 
that  in  Experiments  A  and  B. 

4.2.  Results 

In  all  three  experiments,  mountain  waves  entered  the  mid¬ 
dle  atmosphere  after  a  few  hours  and  produced  overturning 
isentropes  after  ~6  hours.  We  focus  on  temperatures  since 
CRISTA  measured  temperature  perturbations.  Plate  2  plots 
potential  temperatures  0(x,{,t)  and  temperature  perturba¬ 
tions  after/  =  18  hours  over  the  orographic  region 

in  each  experiment.  The  critical  levels  in  Experiments  B  and 
C  efficiently  absorb  mountain  waves,  and  yield  rapid  wave 
shortening  and  intense  breaking  at  closely  underlying  alti¬ 
tudes.  Despite  the  different  wind  profiles  and  the  significant 
amounts  of  overturning  and  mixing  by  this  time,  the  large- 
scale  wave  structures  below  35  km  are  quite  similar  in  each 
experiment. 

The  wave  fields  evolve  significantly  with  time,  as  illus¬ 
trated  in  Figure  2  using  results  from  Experiment  C.  The 
middle  atmosphere  is  dominated  initially  by  short  horizontal 
wavelengths,  which  trigger  vigorous  overturning  and  mixing 
by  8  hours.  The  turbulent  zones  persist  and  seem  to  move 
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Plate  1.  Topographic  elevations  over  southern  South  America,  plotted  as  (a)  filled  contours  in  Cartesian  coordinates, 
(b)  three-dimensional  elevations  in  Mercator  coordinates.  Squares  labeled  1,2  and  3  in  (a)  show  locations  of  temper¬ 
ature  profiles  acquired  by  CRISTA.  Raw  elevations  h(x)  along  the  perpendicular  axis  in  (a)  are  plotted  in  (b).  The 
blue  dot  in  both  figures  is  the  upstream  point.  See  text  for  further  details. 
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Plate  2.  Temperature  fields  after  18  hours  in  Experiments  A,  B  and  C.  Contours  show  isentropes  0(x,  C,  <)•  with 
constant  logarithmic  separation  between  adjacent  contours.  Temperature  perturbations  T'(x,  C,  f )  are  overlayed  using 
the  blue-red  color  scheme  shown  to  the  right.  Underlying  topography  is  shown  in  black. 
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Experiment  C:  t  =  6  hours  Experiment  C:  t  =  8  hours  Experiment  C:  t  =  1 0  hours  Experiment  C:  t  =  15  hours  Experiment  C:  t  =  24  hours 
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Figure  2.  Isentropes  9(x,  C,  t)  from  Experiment  C  after  t  =  6,  8,  10.  15  and  24  hours.  Same  constant  logarithmic 
separation  between  contours  is  used  in  each  panel.  Underlying  topography  is  shaded. 


downstream  after  10  hours.  Thereafter,  progessively  longer 
horizontal  wavelengths  come  to  dominate  the  wave  field,  and 
less  vigorous  overturning  is  evident.  After  24  hours,  waves 
at  lower  stratospheric  levels  are  suppressed  compared  to  their 
earlier  intensities. 


5.  ANALYSIS 
5.1.  Wavelength  Selection 

Preferred  vertical  wavelengths  are  evident  in  Plate  2  and 
Figure  2.  The  vertical  wavelength  of  a  stationary  hydrostatic 
mountain  wave  is  given  theoretically  by 


5.2.  Time  Evolution 

Figures  2  and  3  show  (hat  the  dominant  horizontal  wave¬ 
lengths  change  with  time.  Short  Ar  values  dominate  the  mid¬ 
dle  atmosphere  initially,  giving  way  to  progressively  longer 
horizontal  wavelengths  at  later  times.  The  vertical  group 
velocity  of  any  hydrostatic  wave  (A? /\\  <  1)  is  given  by 


_  X2zN(z) 

91  dm  2irXx 


(4) 


From  (3),  all  waves  below  35  km  have  similar  vertical  wave¬ 
lengths  A.,  and  so  the  vertical  group  velocities  of  various 


X,(z) 


2ttu(z) 

w 


(3) 


where  N  is  Brunt-Vaisala  frequency  and  t  has  its  usual 
meaning  here.  At  15-30  km,  ti(z)  =  20  m  s-1  and 
N  =  0.021  rad  s-1,  which  yields  X2  =  6  km,  in  agreement 
with  the  dominant  vertical  wavelengths  found  below  35  km 
in  all  three  experiments.  In  Experiments  B  and  C,  where 
westward  shear  was  introduced  above  35  km.  A,  decreases 
with  height  in  accordance  with  (3). 

Preferred  horizontal  wavelengths  Xx  also  emerged  in  each 
experiment  (Figure  2).  We  studied  this  spectrally  by  re- 
gridding  T'(*.  C ,  t)  fields  onto  a  regular  height  grid  z,  then 
computing  two-dimensional  power  spectra  Ft(K ,  M)  over 
the  range  x  =  250-1250  km  and  z  =  12-32  km.  As  ex¬ 
pected  from  (3),  for  all  A'  a  sharp  spectral  peak  occurred 
at  a  vertical  wavenumber  |A/|ss  2tr(6  km)-1,  enabling  us  to 
study  overall  spectral  variability  in  these  simulations  using 
the  one-dimensional  spectral  slice  Fj(K,  M  as  A/). 

Spectra  from  each  experiment  are  plotted  in  Figure  3  after 
6,  15  and  24  hours.  The  power  spectrum  of  the  topogra¬ 
phy  h(x)  is  plotted  in  Figure  3d,  with  evident  peaks  marked 
by  light  solid  lines.  These  same  lines  are  also  plotted  in  the 
other  panels  of  Figure  3.  and  show  that  peaks  in  stratospheric 
temperature  spectra  correspond  quite  well  with  peaks  in  the 
topographic  spectrum.  Thus,  horizontal  wavelength  selec¬ 
tion  seems  to  be  governed  by  spectral  features  in  the  Andean 
topography. 
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Figure  3.  One-dimensional  sections  of  the  two-dimensional  power 
spectra  Ft(K,  M)  at  A/  =  Af  =  1/6  eye  km-1  from  Experiments 
A.  B  and  C  after  (a)  6  hours,  (b)  15  hours,  and  (c)  24  hours. 
Spectra  were  computed  from  T‘(x,  z,  t)  fields  in  the  region  z  =  12- 
32  km  and  z  =250-1250  km.  Since  A/  straddles  the  third  and 
fourth  harmonics  of  these  spectra,  we  plot  averages  of  the  slices  at 
these  two  harmonics.  Panel  (d)  plots  power  spectral  density  of  the 
topographic  elevation  h(z).  Local  maxima  in  (d)  are  highlighted 
with  gray  vertical  lines,  and  are  also  shown  in  the  preceding  panels. 
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mountain  wave  modes  are  inversely  proportional  to  their  hor¬ 
izontal  wavelengths.  Ax. 

Using  (3)  and  (4),  a  wave  of  A*  =  50  km  has  a  group 
velocity  «  2.5  m  s-1,  and  thus  takes  ~3.5  hours  to 
propagate  to  z  =  30  km.  Spectra  in  Figure  3a  confirm 
that  the  peak  near  Ax  =  50  km  is  fully  developed  after  6 
hours.  Conversely,  a  Ax  =  150  km  wave  takes  — 10  hours  to 
propagate  to  30  km.  Stratospheric  mountain  wave  energy  at 
this  scale  does  not  develop  fully  until  ~10-15  hours  into  the 
simulation  (Figure  3b).  Note  that  u(z)  was  ramped  up  from 
zero  during  the  first  two  hours,  and  was  <20  m  s-1  below 
17  km,  making  actual  propagation  times  somewhat  longer. 

The  mean  width  of  the  topography  h(x)  is  ~ 300-500  km 
(see  Plate  1 ),  which  might  be  expected  to  force  a  Ax  ~  300- 
500  km  wave  most  strongly.  Indeed,  ridge  widths  here  are 
similar  to  those  in  northern  Scandinavia,  where  very  long 
mountain  waves  often  dominate  the  stratospheric  wave  field 
[Dornbrack  el  al.,  1998,  1999].  Using  (3)  and  (4),  after 
24  hours  a  400  km  wave  would  have  propagated  to  a  height 
of  only  20  km  or  so,  and  thus  would  not  be  fully  developed  in 
the  middle  atmosphere.  Spectra  after  24  hours  confirm  this: 
while  Experiments  A  and  B  show  a  peak  at  a  wavelength 
~400-500  km  in  Figure  3c,  spectra  from  Experiment  C  peak 
most  strongly  nearer  200-250  km.  Thus,  longer  simulations 
are  needed  to  assess  whether  a  Ax  ~400  km  wave  eventually 
dominates  the  middle  atmospheric  wave  fluctuations  [see 
also  Dornbrack  el  al.,  1999], 

Once  each  spectral  peak  has  developed  in  Figure  3,  its 
intensity  tends  to  wane  thereafter,  due  to  wave  breaking  and 
increasing  turbulent  damping.  As  longer  waves  subsequently 
enter  the  region  and  break,  the  additional  turbulence  and 
interactions  they  engender  probably  act  to  damp  these  shorter 
horizontal  wavelengths  further. 

5.3.  Downward  Feedback 

Although  Experiments  A  and  B  used  identical  background 
atmospheres  below  35  km,  the  results  in  Plates  2a  and  2b 
after  18  hours  are  identical  only  below  ~20  km,  and  show 
significant  differences  above  ~25  km.  The  time  evolution 
of  these  differences  can  be  gleaned  from  Figure  3.  After 
6  hours  the  spectra  in  Experiments  A  and  B  are  identical,  but 
show  differences  after  1 5  hours.  By  24  hours  (Figure  3c),  the 
differences  are  as  large  in  places  as  those  from  Experiment  C. 

These  differences  become  evident  below  35  km  after 
waves  start  breaking.  Implicit  numerical  diffusion  in  this 
model  is  very  small,  so  these  differences  must  originate  from 
downward  feedback  of  differing  dynamics  originating  above 
35  km,  where  different  u(z)  profiles  yield  different  wave  and 
turbulence  fields.  Turbulence  itself  cannot  be  advected  this 
far  downwards,  since  we  include  no  mean  vertical  wind. 

Two-dimensional  spectra  Ft(  A',  M)  separate  upward  and 
downward  propagating  waves  [e.g.,  Bacmeister  and  Schoe- 
berl,  1989]:  only  the  upward  components  are  plotted  in  Fig¬ 
ure  3.  After  12-24  hours,  significant  downward-propagating 


spectral  density  appeared  in  each  experiment,  most  promi¬ 
nently  at  long  horizontal  wavelengths.  Vertical  velocity  spec¬ 
tra  Ft(K<  M)  also  showed  downward  propagating  variance 
at  shorter  horizontal  scales.  The  appearance  of  these  features 
in  Experiment  A  indicates  that  this  is  not  a  critical-level  re¬ 
flection  effect,  as  expected  since  mean  Richardson  numbers 
are  too  large  [Jones,  1 968].  Rather,  these  findings  are  in  qual¬ 
itative  agreement  with  previous  simulations  of  secondary 
wave  generation  by  mountain  wave  breaking  [Bacmeister 
and  Schoeberl,  1989 ,  Satomura  and  Sato,  1999], 

Bacmeister  and  Schoeberl  [1989]  found  that  their  down¬ 
ward  propagating  secondary  waves  eventually  led  to  reduc¬ 
tions  in  overall  wave  intensity  well  below  the  regions  of  wave 
breaking.  Attenuated  wave  activity  at  lower  levels  also  arose 
in  our  simulations  (Figure  2e).  Whether  the  attenuation  noted 
here  results  mostly  from  the  effects  of  downward  propagating 
secondary  waves,  or  from  breakdown  of  and  additional  tur¬ 
bulent  damping  by  progressively  longer  upward-propagating 
primary  waves,  is  unclear  at  present. 

5.4.  Comparison  with  CRISTA  Profiles 

The  general  form  of  the  profiles  over  the  Andes  compares 
favorably  with  the  CRISTA  data.  Figure  4  plots  a  sam¬ 
ple  temperature  perturbation  profile  from  Experiment  C  at 
x  =  845  km,  a  location  between  points  1  and  3  in  Plate 
la.  The  CRISTA  profiles  at  these  points  are  overlayed  [see 
Plate  4  of  Gmssmann,  2000].  The  wavelength,  amplitude 
and  general  height  variation  of  the  model  profile  are  all  quite 
similar  to  the  data.  Oscillations  in  Experiment  C  seem  to 
penetrate  slightly  higher  prior  to  shortening  and  dissipating 
below  the  critical  layer.  Note  too  that  the  amplitudes  of  the 
CRISTA  profiles  are  probably  underestimated  due  to  instru¬ 
mental  effects,  and  that  CRISTA  cannot  resolve  any  waves 
of  A.  4:  3-5  km  [Eckermann  and  Preusse,  1999]. 

Model  profiles  from  locations  directly  over  points  1  and 
3  are  not  in  phase  with  the  CRISTA  perturbations.  This 
is  due  to  somewhat  different  horizontal  wavelengths  in  the 
model  simulations.  Eckermann  and  Preusse  [1999]  showed 
that  the  oscillations  measured  by  CRISTA  had  a  horizontal 
wavelength  of  either  400  km  or  400/3  ss  130  km.  Peaks 
somewhat  near  these  values  appear  in  Figure  3c,  although 
other  peaks  also  occur  and  there  are  variations  among  the 
experiments. 

6.  SUMMARY  AND  DISCUSSION 

We  have  developed  a  high-resolution  high-altitude  nonlin¬ 
ear  numerical  model  of  compressible  atmospheric  flow  over 
arbitrary  topography,  and  have  used  it  here  to  simulate  moun¬ 
tain  waves  in  the  middle  atmosphere  over  the  southern  Andes. 
The  results  have  confirmed  that  stratospheric  temperature 
perturbations  measured  by  CRISTA  on  6  November,  1994 
were  produced  by  intense  breaking  mountain  waves  forced 
by  flow  over  the  southern  Andes,  as  argued  by  Eckermann 
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Figure  4.  Black  curve  shows  a  vertical  profile  of  temperature 
perturbations  at  x  =  845  km,  15  hours  into  Experiment  C.  Two 
shaded  solid  curves  reproduce  the  temperature  perturbation  pro¬ 
files  acquired  by  CRISTA  at  locations  1  and  3  in  Plate  la  (after 
Eckermann  and  Preusse  [  1 999]). 

and  Preusse  [1999].  Amplitudes  and  vertical  wavelengths 
from  CRISTA  compare  fairly  well  with  model  profiles.  The 
results  showed  that  attenuated  oscillations  above  30  km  were 
the  result  of  vertical  wavelength  reductions  and  intense  wave 
breaking,  as  conjectured  by  Eckermann  and  Preusse  [1999], 
Recent  analysis  of  satellite  data  by  McLandress  et  al.  [2000] 
has  shown  that  temperature  fluctuations  appear  recurrently  in 
the  middle  atmosphere  over  southern  South  America.  Thus 
these  breaking  mountain  waves  may  be  relevant  to  the  clima¬ 
tology  of  this  region  of  the  middle  atmosphere  [Garcia  and 
Boville,  1994]. 

Presently  we  have  a  limited  understanding  of  the  ampli¬ 
tudes  and  wavelengths  of  mountain  waves  radiated  into  the 
middle  atmosphere  by  complex  topography.  Existing  pa- 
rameterizations  make  use  of  topographic  elevation  spectra, 
focusing  either  on  dominant  peaks  [e.g.,  Bacmeister ,  1993] 
or  their  power-law  shapes  [e.g.,  Shutts,  1995]  to  specify  ra¬ 
diated  mountain  wave  fields.  Our  results,  though  limited  to 
a  specific  two-dimensional  section  of  the  Andes,  produced 
stratospheric  wave  fluctuations  with  distinct  vertical  and  hor¬ 
izontal  wavelengths.  These  horizontal  wavelengths  were 
closely  related  to  peaks  in  the  power  spectrum  of  the  topo¬ 
graphic  elevations  h(x),  lending  some  support  to  spectral 
parameterization  methods.  However,  the  various  horizontal 
wavelengths  evolved  differently  with  time  in  our  model,  due 
mostly  to  varying  vertical  group  velocities.  Wave  fields  were 
further  complicated  at  later  times  by  breaking  and  turbulence, 
leading  not  only  to  increased  damping  of  primary  waves  but 
apparent  generation  of  secondary  waves.  Even  after  24  hours, 
when  simulations  were  terminated  due  to  computational  con¬ 
straints,  mountain  waves  in  the  middle  atmosphere  were  still 
evolving. 

Actual  topography  and  flow  fields  here  are  intrinsically 
three-dimensional  [Plate  lb;  Figure  la;  Seluchi  etal.,  1998], 
which  can  influence  mountain  waves  significantly  [Shutts, 


1998].  While  two-dimensional  models  with  turbulence  pa- 
rameterizations  provide  reasonable  simulations  of  breaking 
gravity  waves  [Liuetal.,  1999;  Doyle  etal.,  2000],  long-term 
fine-scale  simulations  of  wave  breakdown  beneath  a  critical 
level  seem  to  require  three-dimensional  models  [Winters and 
D'Asaro,  1994;  Dombrack ,  1998;  Fritts  and  Weme,  2000], 
Additionally,  the  model  here  did  not  include  mean  flow  evo¬ 
lution,  which  may  be  important  in  this  problem.  Thus,  multi¬ 
day  interactive  runs  using  three-dimensional  models  may  be 
needed  to  describe  these  evolving  mountain  wave  fields  fully. 
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